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ABSTRACT 



Aims. We perform a systematic combined photometric and kinematic analysis of a sample of globular clusters under different re- 
laxation conditions, based on their core relaxation time (as listed in available catalogs), by means of two well-known families of 
spherical stellar dynamical models. Systems characterized by shorter relaxation time scales are expected to be better described by 
isotropic King models, while less relaxed systems might be interpreted by means of non-truncated, radially-biased anisotropic / <v) 
models, originally designed to represent stellar systems produced by a violent relaxation formation process and applied here for the 
first time to the study of globular clusters. 

Methods. The comparison between dynamical models and observations is performed by fitting simultaneously surface brightness and 
velocity dispersion profiles. For each globular cluster, the best-fit model in each family is identified, along with a full error analysis 
on the relevant parameters. Detailed structural properties and mass-to-light ratios are also explicitly derived. 

Results. We find that King models usually offer a good representation of the observed photometric profiles, but often lead to less 
satisfactory fits to the kinematic profiles, independently of the relaxation condition of the systems. For some less relaxed clusters, 
/ (v) models provide a good description of both observed profiles. Some derived structural characteristics, such as the total mass or 
the half-mass radius, turn out to be significantly model-dependent. The analysis confirms that, to answer some important dynamical 
questions that bear on the formation and evolution of globular clusters, it would be highly desirable to acquire larger numbers of 
accurate kinematic data-points, well distributed over the cluster field. 

Key words. (Galaxy:)globular clusters: general - (Galaxy:)globular clusters:individual: NGC 104 (47 Tuc), NGC 288, NGC 362, 
NGC 2419, NGC 3201, NGC 5139 (co Cen), NGC 6121 (M4), NGC 6218 (M12), NGC 6254 (M10), NGC 6341 (M92), NGC 6656 
(M22), NGC 6809 (M55), NGC 7078 (M15) 



1. Introduction 

It is well known that for most globular clusters the relevant relax- 
ation times are shorter than their age. Therefore, since the two- 
body relaxation processes have had time to act, these systems 
may be considered to be close to thermodynamical relaxation, 
with a distribution function that should be close to a Maxwellian. 



The spherical King models (King 1966) are defined by a 
distribution function of this type, with a truncation in energy 
introduced to mimic the role of tides (for the relevant defi- 
nition and a summary of the main properties, see App. IA. lb . 
In general, they are thought to be a physically justified sim- 
ple tool to represent quasi-relaxed stellar systems, in partic- 
ular well suited to describe globular clusters. In this context, 
King models are usually considered as the correct zeroth-order 
dynamical r eference model. The available catalogs (e.g., see 
Harrisl l20Tol) record a number of structural properties of glob- 
ular clusters that are derived from an interpretation of observa- 
tional data based on the use of King models. The family of King 
models is also frequently used in dynamical studies aimed at 
identifying given physical mechanisms and represents a most 
popular choice of initial conditions for numerical simulations 
designed to in vestigate the dynamical ev olution of star clus- 



ters (e.g., see IChernoff & Weinberg! [l 990. and, more recently, 



iGieles & Baumgardtll2008l) 

In reality, this simple physical picture suffers from a number 
of limitations, which become more evident now that much im- 



proved observations allow us to go well beyond a zeroth-order 
picture of globular clusters. Of course, numerical simulations 
can now model in a realistic way the important effects of dis- 
creteness, mass segregation, bin aries, rotation, core c ollapse, and 
time-dependent tides (e.g., see iKiipper et al.|[2010h . which are 
known to govern the evolution of globular clusters. Still, in rela- 
tion to the general issue of relaxation it would be interesting to 
test further to what extent King models (or alternative idealized 
models) actually work as a useful paradigm to represent small, 
quasi-relaxed stellar systems. 

From the observational point of view, there remains indeed 
one unsatisfactory point, which this paper would like to ad- 
dress, at least on the basis of the currently available data. This 
is that the application (and thus the test) of the spherical King 
models is frequently carried out, in pra ctice, only in rel ation 
to the available photometric profile s (see lTrager et alJ[l995l and 
iMcLaughlin & van der M arel 2005), without the corresponding 
tests on the associated kinematical profiles, in contrast to what 
is routinely done in the studies of early-type galaxies. Some 
studies of Galactic globular clusters, constrained simultaneously 
by density and velocity dispersion profiles, are actually avail- 
able, mainly based on multi-mass M ichie-Kin g models (e.g., see 
i Gunn & Griffinlll979llMevlan et alj|1995l and lMevlan & Mavorf 
Il99ll for M3. id Centauri, and NGC 6397, respectively), but a 
systematic and homogeneous investigation is still missing. Note 
that, on the galactic side, deep investigations of these issues, 
starting with the mid 1970s, have led to the remarkable discov- 
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ery that bright ellipticals are generally supported by anisotropic 
pressure and contain significant amounts o f dark matter inside 
the effective radius R e (e.g., see Chap. 24 in Bertin 2000). From 
the study of elliptical galaxies it has also been learned that 
structurally different models (as diagnosed by their kinemat- 
ics or characterized by their virial coefficients) may have re- 
markably similar photometric profiles (e.g., see Appendix D of 
iBertin et alj|2002l) . 

In this paper, we wish to check whether indeed King models 
perform better in more relaxed systems, by studying the com- 
bined photometric and kinematical profiles for a sample of glob- 
ular clusters . We will refer to the surface brightness profiles 
collected by iTrager et al.l (1 19951) (supplemented by more recent 
data, when available), in order to deal with a homogeneous sam- 
ple, and limit our discussion to one-component dynamical mod- 
els (leaving aside the issues of core collapse and mass segrega- 
tion; we will thus exclude from our sample clusters with evi- 
dence of core collapse). Our sample will then be basically de- 
fined by the requirement that a sufficient number of kinematical 
data-points is available from the literature, so as to define a rea- 
sonably accurate and radially extended kinematic profile for a 
test of a dynamical model at the global level. 

Many large globular clusters have very long relaxation times. 
Therefore, we have decided to mode l the same data als o by 
means of models (the models; see lBertin & Trent! I2003h ex- 
plicitly constructed for the context of violently relaxed elliptical 
galaxies. In other words, we wish to check whether less relaxed 
clusters tend to conform to the picture of formation via incom- 
plete violent relaxation, which has the characteristic signature 
of radially-biased pressure anisotropy for less bound stars. The 
use of the / <v) models is preferred to o ther options ( for exam- 
ple to the use of King-Michie models, Michie 1963), because 
these models are based on a detailed physical justification and 
have been shown to perform well both in relation to the obser- 
vations of bright ellipticals and to the properties of the products 
of incomplete violent relaxation found in numerical simulations 
of collisionless collapse (over a range of nine orders of mag- 
nitude in the computed density profiles, with an excellent fit to 
the properties of the generated pressure anisotropy profiles; see 
Trenti, Bertin, van Albada 2005). For a fair comparison with the 
King models we should have referred to truncated / <v) models; 
for the present simple exploratory investigation, and in order to 
keep the comparison between models characterized by the same 
number of parameters, we have decided to use the non-truncated 
f {v) models. 

The paper is organized as follows. In Sect.[2]we introduce the 
sample of globular clusters selected for this study. In Sect. [3] we 
describe the available data sets and the procedure we followed 
to obtain the radial profiles used in our analysis. In Sect. [4] we 
show the results of our work and in Sect. [5] we draw our conclu- 
sions. A brief description of the characteristics of the dynamical 
models is provided in Appendix lAl in Appendix[B]we describe 
the adopted fitting procedure. 

2. The selected sample 

In this paper we wish to analyze globular clusters characterized 
by different relaxation conditions, measured by the central (core ) 
relaxation time T c as listed in the Harris catalog (Harris 2010). 
We decided to order globular clusters according to the core re- 
laxation time, rather than to the half-mass relaxation time, be- 
cause in general T c is less model-dependent (as pointed out in 
Sect. 14.21 ). We checked that, by ordering the sample with respect 
to the half-mass relaxation time, some changes in the composi- 



tion of the three classes identified below would occur (this fact 
can be easily seen by inspecting the values of these time scales 
for the globular clusters listed in TablefTJ, but with minor effects 
in relation to the general conclusions of the paper. 

Looking at the values of the relevant relaxation times 
recorded in the available catalog for the globular cluster system 
of our Galaxy, we see that the distinction between relaxed and 
partially relaxed globular clusters is not sharp. Therefore, we in- 
troduce a simple criterion to classify globular clusters according 
to their relaxation state. We consider three relaxation classes: 
relaxed globular clusters, for which log T c < 8 (first class); glob- 
ular clusters in an intermediate relaxation condition, for which 
8 < log T c < 9 (second class); partially relaxed globular clus- 
ters, for which log T c > 9 (third class; T c is expressed in years). 

The sample of globular clusters has been selected on the ba- 
sis of the following criteria: (i) We exclude post-core-collapse 
globular clusters, that is , we reject the clusters labeled as post- 
core-collapse by lHarrisI d2010h . The reason is that we wish to 
test King models on the global scale, avoiding the subtle mod- 
eling issues that characterize the central regions of these sys- 
tems, especially if phenomena a priori known to go beyond the 
King modeling are involved, (ii) We choose globular clusters 
for which an accurate and extended surface brightness profile 
is available, (iii) We select clusters for which at least 140 stel- 
lar radial (line-of-sight) velocities have already been measured. 
We impose a lower limit to the number of measured velocities, 
because we wish to extract from such data a reasonably well- 
defined velocity dispersion profile; the value of this limit is fixed 
in a way that allows us to include in the analysis globular clusters 
that belong to the different relaxation classes defined previously, 
(iv) We exclude from our list all the clusters that have less than 
35 stellar radial velocities inside the projected half-light radius 
R e (also called effective radius; the values of this quantity are 
reported in the Harris catalog, where the notation is rh). We in- 
troduce this further requirement, because we wish to analyze ve- 
locity dispersion profiles that characterize the stellar systems on 
the largest radial extent. For bright elliptical galaxies the kine- 
matical data-points inside R e are usually the easiest to get, and 
often turn out to discriminate among very different models. To 
test how well the King models perform, the central regions are a 
natural ground for comparison with other models. 

The most restrictive elements in identifying a significant 
sample of globular clusters are the requirements on the radial 
velocity data, because for only few globular clusters the desired 
data are available. Indeed, of the 28 Galactic globular clusters 
with a reasonable number of radial velocities (i.e., at least 40 
line-of-sight velocity measures on the entire spatial extent of the 
cluster), 3 are flagged as post-core-collapse, 9 have less than 140 
velocity data, and 3 have less than 35 data inside their projected 
half-light radius. In the Harris catalog, NGC 362 and NGC 7078 
are indicated as possible post-core-collapse clusters, but we de- 
cided to keep them because, according to their concentration pa- 
rameter in the Harris catalog, it is still possible to obtain an ac- 
ceptable fit with King models. In this way, we are left with 13 
globular clusters that match our selection criteria. In relation to 
the relaxation classes defined above, our set of globular clusters 
contains 5 well-relaxed clusters, 6 clusters in an intermediate re- 
laxation condition, and 2 partially relaxed clusters. 

To better characterize our sample in terms of the radial extent 
of their radial velocity data, we consider the ratio of the radius 
of the last kinematical point to the projected half-light radiufl 



1 We always indicate with R the projected (two-dimensional) radial 
scales, and with r the intrinsic (three-dimensional) radial scales. 
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Table 1. The selected globular clusters. 



Globular Cluster 


d Q 


C 


lo g r c 


logT M 


e 


Vrot 




N e 


1> V 

A K 

R« 


n K 
/1r 


Ref. 


NGC 362 


8.6 


1.76 


7.76 


8.93 


0.01 


0.0 


208 


92 


4.22 


0.33 


(1) 


NGC 7078 (M15) 


10.4 


2.29 


7.84 


9.32 


0.05 


1.7 


Mil 


1298 


16.94 


0.62 


(2) 


NGC 104 (47 Tuc) 


4.5 


2.07 


7.84 


9.55 


0.09 


2.2 


2638 


709 


19.27 


1.44 


(3), (4) 


NGC 6121 (M4) 


2.2 


1.65 


7.90 


8.93 


0.00 


0.9 


200 


55 


10.36 


0.87 


(4) 


NGC 6341 (M92) 


8.3 


1.68 


7.96 


9.02 


0.10 


2.5 


295 


42 


13.96 


1.14 


(5) 


NGC 6218 (M12) 


4.8 


1.34 


8.19 


8.87 


0.04 


0.15 


242 


58 


10.38 


1.06 


(4) 


NGC 6254 (M10) 


4.4 


1.38 


8.21 


8.90 


0.00 




147 


47 


5.22 


0.55 


(6) 


NGC 6656 (M22) 


3.2 


1.38 


8.53 


9.23 


0.14 


1.5 


345 


116 


8.40 


0.88 


(4) 


NGC 3201 


4.9 


1.29 


8.61 


9.27 


0.12 


1.2 


399 


201 


10.35 


1.27 


(7) 


NGC 6809 (M55) 


5.4 


0.93 


8.90 


9.29 


0.02 


0.25 


728 


311 


7.79 


1.44 


(4) 


NGC 288 


8.9 


0.99 


8.99 


9.32 




0.25 


171 


68 


5.53 


0.93 


(4), (6) 


NGC 5139 (w Cen) 


5.2 


1.31 


9.60 


10.09 


0.17 


7.9 


2060 


554 


5.97 


0.62 


(8), (9) 


NGC 2419 


82.6 


1.37 


9.87 


10.63 


0.03 


0.6 


166 


38 


14.63 


1.74 


(10) 



Notes. From left to right, the following quantities are displayed: distance from the Sun (kpc), concentration parameter, logarithm of the core 
relaxation time (years), logarithm of the half-mass relaxation time (years), ellipticity, rotational velocity (km s~'), total number of velocity data- 
points available, number of velocity data-points inside the projected half-light radius, ratio of the radius of the outermost velocity point to the 
projected half-light radius, and ratio of the radius of the outermost velocity point to the truncation radius. The sources of the kinematical data are 
listed in the last column (see main text for references of the other quantities). 

References. (1) iFischer et alJ 1 19931: (2 ) IGebhardt et all l2000l : (3) IGebhardt et all 119951: (4) lLaneetal.1 1201 It (5) iDrukier et all 120071: (6) 
iCarretta et a l. 2009; (7) ICote et al .11 19951: (8) lMavor et alJ 19971 : (9) lReiins et alj2006k (lO) IIbataet alj|2011l . 



Rk/Rs, and the ratio of the radius of the last kinematical point to 
the truncation radiufl R^/r^. We judge the following values of 
the two ratios, Rk/Ri: ^ 3 and R&lr a > 0.8, to be satisfactory. 
All the selected globular clusters satisfy the first relation, and all 
but four globular clusters satisfy the second condition. 

Table [TJgives the sample of selected globular clusters, listed 
in order of increasing core relaxation time log T c . The first part 
of the table contains relaxed globular clusters, the second part 
those in an intermediate relaxation condition, and partially re- 
laxed clusters are shown in the last part. For each object, the val- 
ues of the adopted cluster distance from the Sun cIq (expressed 
in kpc), the concentration parameter C, the logarithm of the core 
relaxation time log T c , the logarithm of the half-mass relaxation 
time log Tm (where T c and Tm, in other papers often indicated 
with the symbols f rc and f r h, are expressed in years) and the el- 
lipticity e are recorded (as listed in the Harris 2010 catalog). In 
addition, the maximum rotati onal velocity v rot (in km s _1 ; the 
references for these values are lLane et aTll201 ll for the globular 
clusters f or which we use the kin ematic data published in this 
paper, and Me vlan & Heggid 19971 for the others), the number of 
velocity data-points available N v , the number of velocity data- 
points inside the projected half-light radius N e , the ratios of the 
radius of the outermost velocity point R& to the projected half- 
light radius R e (in the Harris catalog and in other papers often 
indicated as r^) and to the truncation radius r tr , and the sources 
of the kinematic data are given in the last columns. 

For completeness, in Table|2]we list globular clusters with a 
fairly rich kinematical information available (line-of-sight veloc- 
ity measurements for at least 40 stars), but rejected in the present 
analysis because they do not satisfy our selection criteria. The 
table is organized in three parts, thus showing which criterion is 
not met: in the first part we list post-core-collapse globular clus- 

2 In the latest version of the Harris catalog (Harris 2010|), the val- 
ues of the truncation radii are not listed; we calculated them from the 
available values of R c and C (see Appendix I A. lb , as: r tr = R c 10 c (see 
notes in the Harris catalog bibliography). For the three core-collapsed 
globular clusters in Table [2] for which the concentration parameter is 
arbitrarily set to C = 2.50, the value of r tr should be considered only as 
indicative. 



ters, in the second part those with an insufficient total number 
of kinematic data {N w ), and in the third part those with an insuf- 
ficient number of data inside the half-light radius (N e ). In each 
part of the table, globular clusters are listed in order of increas- 
ing number of total stellar velocity data N v . Column entries are 
the same as in Table[T](except for cluster distance, ellipticity, and 
rotational velocity, which are not recorded). 



3. The data sets 

3.1. The surface brightness profiles 

To deal with a homogeneous sample , we decided to use the sur- 
face brightness profiles pr ovided bvlTrager et alJ (I1995I) (this i s 
the same starting point of iMcLaughlin & van der Marell l2005). 
This choice guarantees that the profiles have been constructed 
with the same method, even though the actual data come from 
different sources. 

For each globular cluster, the profile is composed of A^ p pho- 
tometric data-points, given by log /?,-, the logarithm of the radius 
Ri, measured in arcsec, and by my(Ri), the V band surface bright- 
ness measured in mag arcsec -2 at the radial position /?,-. The data 
also include myc(Ri), the surface brightness calculated with the 
Chebyshev polynomials, which provides an accurate approxi- 
mation to the overall profile; my(Rj) - myc(Ri), the Chebyshev 
residual; w,-, a weight that the authors assign to each measure- 
ment. All the surface brightness profiles considered in this paper 
are relative to the V band. 

These data have to be properly corrected and treated before 
a comparison can be made with the theoretical models. First, 
we introduced an extinction correction, under the assumption 
that such extinction can be considered to be constant over the 
entire extent of each globular cluster. We calculated the extinc- 
tion Ay from the reddening listed in the Harris catalog and ob- 
tained m(Ri) = my{Rj) - Ay, for i = 1, . . .,N V . This is the 
only c orrection applied to the data: we assume that lTrager et alJ 
(1 19951) already removed any foreground and background con- 
tamination that could affect the measurements. Then, we fol- 
lowed the procedure described by McLaugh lin & van der Mare! 
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Table 2. Rejected globular clusters. For a description of column entries, see TableQ] 



Globular Cluster 


C 


io g r c 


logT M 




A'e 


7T 


ro- 


Ref. 


NGC 6397 


2.50 


4.94 


8.60 


103 


103 


0.37 


om 


(1) 


NGC 7099 (M30) 


2.50 


6.37 


8.88 


196 


33 


11.20 


0.61 


(2) 


NGC 6752 


2.50 


6.88 


8.87 


437 


75 


15.24 


0.54 


(2) 


IC 4499 


1.21 


9.21 


9.73 


43 


11 


6.93 


0.87 


(3) 


NGC 6171 (M107) 


1.53 


8.06 


9.00 


66 


37 


3.67 


0.40 


(4) 


NGC 5466 


1.04 


9.35 


9.76 


66 


49 


1.69 


0.20 


(5) 


NGC 7089 (M2) 


1.59 


8.48 


9.40 


69 


19 


8.02 


0.68 


(6) 


NGC 5053 


0.74 


9.81 


9.87 


71 


27 


3.28 


0.75 


(7) 


NGC 4590 (M68) 


1.41 


8.45 


9.27 


122 


29 


6.69 


0.68 


(8) 


NGC 2808 


1.56 


8.24 


9.15 


123 


7 


8.41 


0.74 


(9) 


NGC 6205 (M13) 


1.53 


8.51 


9.30 


123 


10 


7.92 


0.64 


(10) 


NGC 5904 (M5) 


1.73 


8.28 


9.41 


136 


18 


6.61 


0.50 


(8) 


NGC 1904 (M79) 


1.70 


7.83 


8.95 


146 


2 


14.46 


1.17 


(11) 


NGC 5024 (M53) 


1.72 


8.73 


9.76 


180 


4 


16.24 


1.16 


(2) 


NGC 1851 


1.86 


7.43 


8.82 


184 





21.21 


1.66 


(11) 



References. (1) Geb hardt et al.ll995k (2)lLane et all20Ul; (3)lHankev & Cold20To|; (4)lPiatek et aljl9 94: (5) Sh etrone et alj201(t (6) lPrvor et all 
119861 : (7) 1 Yan & Cohenll 19961 : (8) ICarretta et al.ll2009b (9) ICarretta et al.ll2006l : (lO) lMeszaros et alJ2009l: (ll) IScarpa et alj201lL 



(2005) to estimate the uncertainties &m t on the data: starting from 
the weights, we calculatecQ dm, — cr^/wi (note that in the fits we 
use only the poi nts with weights w, > 0. 15 in the orig inal profile, 
as suggested by McLaug hlin & van der Marel 2005b. 

Usually, t he les s reliable parts in the profiles published by 
iTrager et all (1 19951) are the central regions. Therefore, we de- 
cided to combine these profil e with the more recen t and a ccurate 
surface brightness profiles bv lNovola & Gebhardtl (T20 06). when 
available. For the globular clusters NGC 104 and NGC 6341, 
we simply combined the profiles from the two sources; for NGC 
6254 and NGC 7078, we decided to combine the two data sets by 
removing the points f romlTrager et al.l (1 19951). wh en they do not 
match the profile by iNovola & Gebhardtl (|2006) (the removed 
points are inside 3.5" and 14.5", respectively). In the case of 
NGC 7078, it should be emphasized that, with this treatment, 
the profile changes significantly in the central regions, and the 
central slope becomes steeper. 

In the case of NGC 5139 (w Cen) we decided to add to the 
surface brig htness profile the inn er points that Eva Noyola kindly 
provided us (INovola et al.ll2008l) : still the number of data-points 
in our final composite surface brightness profile of this cluster is 
significantly smaller than that of the other clusters. 



3.2. The velocity dispersion profiles 

We divided the data (see Table[TJfor detailed references) in sev- 
eral radial bins containing an equal number of stars; we chose 
the binning that represents the best compromise between having 
a rich profile and having accurate points (by increasing the num- 
ber of bins, that is by decreasing the number of points per bin, 
the errors on the velocity dispersion increase). In principle, we 
might consider using unbinned data, to a void loss of informatio n 
(for the study of dwarf spheroidals, see Wil kinson et al . 2002). 
Here we preferred to follow the more traditional approach of 
constructing the associated one-dimensional profiles, a method 
that can be applied in a similar way to both kinematical and pho- 
tometric data and that allows us to follow well-established fitting 
procedures used in the past (especially in studies of elliptical 
galaxies). For the clusters with less numerous kinematical data- 



points, we did experiment with changing the radial binning and 
checked the related consequences. 

The method used to calculate the mean velocity and the ve- 
locity dispersio n from stellar radial ve locities is basically the one 
described by Pryor & Meylanj d 1 993b . We started by calculating 
the mean v r and the dispersion <x c for the entire set of velocities. 
The mean velocity represents the overall velocity of the entire 
cluster. Then, taking this value as a constant for the entire clus- 
ter, we calculated the line-of-sight velocity dispersion o-(Rj) and 
the related accuracy £cr, inside the bins in which the data have 
been divided. For each bin, we indicate the distance from the 
center R, as the mean of the radial positions of the stars that it 
contains. In this paper we ignore the possible presence of rota- 
tion and therefore consider the various kinematical data-points in 
each bin, after subtraction of the systemic velocity, to contribute 
only to velocity dispersion (random motions). 

For the majority of the clusters in our sample, only one data- 
set of stellar radial velocities is available; in the following we 
discuss in detail the cases in which a composition of different 
data-sets has been performed or which require some additional 
comments. For NGC 104, two data-sets of radial veloci t ies are 
available; we noticed that the data from Gebha rdt et al. (1 19951) 
are m ore centrally concentrated than those from lLane et al.l 
(201 1)Q In order to have a complete sampling on the entire ra- 
dial extent of the cluster, we decided to defin e a mixed data-set, 
composed of 499 data from Geb hardt et alj ( 19951) . located in- 
side 100", and 2139 data from Lane et alJ (1201 1 ). located out- 
side that radius. The second case is cluster NGC 288, studied by 
ICarretta et al.l d2009l) and by lLane et alj d201 ll) . Since these pa- 
pers publish the coordinates of each star, we were able to single 
out the stars in common between th e two data-sets: for th e stars 
in the overlap, velocity measures bv lCarretta et al.l d2009l) . being 
more accurate, have been preferred. Finally, we excluded three 



3 CTp is a constant that varies from cluster to cluster, the value o f 
which can be found in Table 6 in McLaughli n & van der MarelN2005h . 



4 iGiersz & Heggid J2011I) make some cautionary remarks about the 
velocity dispersion profile repor ted by these au thors. In particular, the 
selection criteria adopted by lLane et al] d201 lh could lead to the ex- 
clusion of some high-velocity stars, with consequent lowering of the 
central velocity dispersion, and to the inclusion of nonmember stars af- 
fecting the outer part of the profile. By using the composite data-set 
described above, we should be able to obtain reliable values of the ve- 
locity dispersion in the central regions, while the outermost points may 
still be affected by the inclusion of nonmember stars. 
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Table 3. Dimensionless parameters and physical scales of the best-fit models. 









King 


Models 












f v) Models 






NGC 






ro 


y"0 




V 














V 


(1) 




(2) 


(3) 


(4) 




(5) 




(6) 


(7) 


(8) 




(9) 


104 


8 


58 + 0.01 


23.09 ± 0.23 


14.33 ± 


0.01 


12.27 ±0.19 


8 


21 


± 0.02 


250.44 ± 1.46 


14.29 ± 


0.01 


14.07 ± 0.22 


288 


4 


82 ±0.10 


91.03+2.86 


20.02 ± 


0.03 


2.85 ±0.19 


3 


91 


± 0.24 


79.54 ± 5.37 


19.88 ± 


0.03 


3.70 ± 0.25 


362 


7 


80 ± 0.03 


10.21 ±0.11 


14.66 ± 


0.01 


8.31 ± 0.44 


6 


86 


±0.05 


52.73 ± 1.09 


14.70 ± 


0.01 


9.26 ± 0.50 


2419 


6 


62 ± 0.04 


19.70 + 0.31 


19.43 ± 


0.03 


5.04 ± 0.43 


4 


24 


±0.09 


29.86 ± 0.83 


19.44 ± 


0.04 


7.28 ± 0.61 


3201 


6 


17 + 0.11 


76.99 ± 3.05 


18.35 ± 


0.08 


4.28 ±0.19 


4 


09 


±0.40 


103.48 ± 11.71 


18.33 ± 


0.08 


5.01 ± 0.23 


5139 


6 


27 ± 0.05 


136.94 ± 2.33 


16.42 ± 


0.04 


14.83 ± 0.25 


4 


31 


±0.07 


150.63 ± 3.07 


16.35 ± 


0.04 


23.41 ± 0.40 


6121 


7 


32 ± 0.07 


74.56 ± 1.76 


17.00 ± 


0.11 


4.01 ± 0.30 


7 


39 


±0.09 


464.49 ± 20.40 


17.01 ± 


0.11 


4.21 ±0.31 


6218 


6 


1 1 ± 0.07 


51.56 ± 1.51 


17.65 ± 


0.07 


3.93 ± 0.30 


4 


00 


± 0.14 


60.69 ± 2.52 


17.57 ± 


0.07 


5.58 ± 0.42 


6254 


6 


26 ± 0.04 


53.63 ± 0.60 


16.88 ± 


0.09 


6.21 ± 0.37 


2 


67 


±0.13 


51.48 ± 1.45 


16.84 ± 


0.09 


9.69 ± 0.59 


6341 


7 


54 ± 0.02 


14.72 ±0.13 


15.31 ± 


0.01 


9.28 ± 0.41 


5 


99 


±0.04 


50.00 ± 0.80 


15.44 ± 


0.01 


12.77 ± 0.56 


6656 


6 


47 + 0.11 


86.18 ±2.34 


16.41 ± 


0.11 


6.47 ± 0.38 


5 


99 


±0.26 


241.58 ± 26.46 


16.41 ± 


0.11 


7.19 ±0.42 


6809 


4 


44 + 0.11 


129.11 ±4.06 


19.12 ± 


0.04 


2.92 ±0.13 


3 


92 


±0.19 


101.53 ±5.39 


18.99 ± 


0.04 


3.61 ±0.17 


7078 


8 


09 ± 0.02 


7.72 ±0.13 


14.07 ± 


0.03 


11.83 ±0.24 


8 


17 


±0.05 


65.88 ± 0.94 


13.59 ± 


0.06 


12.79 ± 0.26 



Notes. For each cluster, named in column (1), for King and / <v) models, we list: the concentration parameter x ¥, in other papers often indicated as 
Wo (Col. (2) and (6)), the scale radius, expressed in arcsec (ro in Col. (3) and r sca i e in Col. (7), as defined in equations dA.2t and JA.5b : note that 
they are intrinsic quantities; they are recorded here in arcseconds, for easier comparison with the observations, as shown in Figs. UI3t , the V band 
central surface brightness /jq in mag arcsec -2 (Col. (4) and (8)), and the central line-of-sight velocity dispersion V in km s _1 (Col. (5) and (9)). 
Formal errors on the various parameters are also recorded (see Appendix lB.4t . 



stars for which the value of the velocity deviates by more than 
4cr c from the mean radial velocity v r , obtaining a final sample of 
171 data. 

For the globular clusters NGC 61 21, NGC 6656 an d NGC 
6809, in addition to recent data from lLane et al. l d20Tlh . older 
radial velocity measures are available in th e literature, but w e 
decided to consider only the data-sets from lLane et all d2011h . 
because they are more complete and more radially extended. 
In the case of NGC 5139, we merged the largest sample of 



velocity data available ( Reiins et al.ll2.Q06b with the sample pro- 
vided by Mavor et al. ( 1997), which provides measurements for 
stars located in the central region of the cluster. A delicate is- 
sue regarding this cluster is the controversial position of its cen- 
ter, which plays an important role also in our analysis, because 
we wish to build a radial-dependent velocity dispersion profile, 
starting from st ellar positions exp ressed in right ascension and 
declination (for lReiins et alJ l2006). To carry out a proper merg- 
ing of the tw o data-sets, we have used the position of the center 
proposed by iMavor et alj d!997|) ; to calculate radial distances, 
we ha ve followed the procedure described by Ivan de Ven et alJ 
(120061) . 



4. Results 

4.1. Relaxation classes 

The values of the dimensionless parameters and the physical 
scales of the two families of models determined by the photo- 
metric and kinematic fits are presented in Table [3] (the fitting 
procedure is described in AppendixlBt. Note that, in general, for 
the / (v) models relatively low values of the concentration param- 
eter *¥ are identified, consistent with the comment given at the 
end of Appendix lA.2l Quantitative information about the best-fit 
models and the observational profiles, such as the number of the 
photometric and kinematic points, the values of the relevant re- 
duced chi-squared, and the corresponding residuals, are listed in 
Table|4] To evaluate the goodness of the fits, in Table[5]we com- 
pare the values of the reduced photometric and kinematic chi- 
squared, denoted by xi and x\ respectively, with the two-sided 
90% confidence interval (CI), calculated with respect to the re- 



duced ^-distribution, characterized by the appropriate number 
of degrees of freedom (see App. lB.3l for details). 

The surface brightness and the line-of-sight velocity disper- 
sion profiles determined by the fit procedure for the models, to- 
gether with the observed profiles for the clusters in the first, sec- 
ond, and third relaxation class, are shown in Figs. [1] El and [3] 
respectively. In the panels, solid lines correspond to the best-fit 
King models and dotted lines to the best-fit / (v) models. The ver- 
tical solid line marks the position of the King model projected 
half-light radius, the dotted one the position of the / (v) model 
projected half -light radius . For t he surface brightness profiles, 
the data from Traa er et alJ d 19951) are indicated with circles, the 
data from other sources with squares. For the velocity disper- 
sion data, the horizontal bars indicate the length of the radial bin 
in which the data-points have been calculated; they do not have 
a role in determining the fit. For each data-point the errors are 
shown as vertical error bars. Note that, even if we insisted on se- 
lecting clusters with a reasonable number of data inside R e , for 
about half of the clusters, the kinematic profiles are undersam- 
pled in their central region. 

In the following part of this subsection we will try to give a 
general assessment of the quality of the fits in the various cases. 
As a general rule, for a given cluster the family of models yield- 
ing the lowest values of the best-fit photometric and kinematic 
chi-squared is preferred. In practice, the objective properties of 
the fits are best obtained by checking directly the values pro- 
vided in Tables 3] and |5] and by inspection of Figs. [T][3] Table [4] 
shows that generally, for a given cluster, x\ > x%, because the 
photometric profiles are characterized by a larger number of 
data-points with reported smaller error-bars. Table[5]shows that 
NGC 6656 is the only cluster for which the values of xi an d 
x\ for the two models are inside the relevant 90% CI. Within 
the King modeling, four clusters have xi inside the 90% CI and 

nine have ~x\ inside the 90% CI. A visual inspection of the fits 
(in particular, see the kinematic fits for NGC 362, NGC 3201, 
NGC 6121, NGC 6656, NGC 6809 and the inner photometric 
profiles of NGC 3201, NGC 6121, NGC 6218, NGC 6656, and 
NGC 7078) also suggests that some systematic trends are clearly 
missed by the two families of models used. 
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Table 4. Quality of the fits. 













King Models 










/W Models 






NGC 








<A^> 


(A/i)max Xt 


(Aa) 


(Acr) max 




<A/i> 


(A/i)n,ax ^ 


<A<r> 


(A<x) max 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) (7) 


(8) 


(9) 


(10) 


(11) 


(12) (13) 


(14) 


(15) 



104 


231 


16 


3.487 


0.41 


4.68 


7.411 


1.03 


3.39 


6.433 


0.24 


1.53 


10.367 


1.34 


2.98 


288 


85 


6 


1.251 


0.30 


0.98 


0.442 


0.21 


0.33 


3.891 


0.46 


1.10 


2.040 


0.46 


0.78 


362 


239 


8 


3.113 


0.58 


7.09 


1.307 


0.99 


1.74 


1.563 


0.15 


0.92 


3.345 


1.55 


2.45 


2419 


137 


6 


1.983 


0.21 


1.10 


1.344 


0.98 


2.21 


1.492 


0.16 


0.83 


0.471 


0.50 


0.83 


3201 


80 


16 


1.308 


0.38 


1.49 


1.783 


0.83 


1.73 


1.289 


0.36 


1.48 


4.005 


1.21 


2.29 


5139 


72 


37 


3.750 


0.36 


2.08 


1.974 


1.73 


4.90 


21.742 


0.80 


1.62 


3.406 


2.07 


4.38 


6121 


228 


10 


1.460 


0.27 


1.33 


0.450 


0.47 


0.94 


1.710 


0.29 


1.33 


0.581 


0.52 


0.93 


6218 


143 


11 


1.185 


0.32 


1.12 


0.584 


0.54 


0.91 


2.663 


0.40 


1.12 


0.765 


0.59 


1.10 


6254 


162 


6 


5.046 


0.37 


2.75 


0.606 


0.48 


0.65 


4.372 


0.22 


1.12 


1.844 


0.89 


1.15 


6341 


118 


8 


8.439 


0.41 


2.51 


1.418 


0.51 


1.02 


20.589 


0.33 


1.06 


2.354 


1.01 


2.37 


6656 


143 


7 


1.019 


0.23 


0.66 


0.942 


0.67 


1.36 


1.056 


0.23 


0.66 


1.699 


0.89 


1.83 


6809 


114 


13 


1.165 


0.32 


1.04 


1.103 


0.40 


0.96 


4.404 


0.59 


1.34 


2.967 


0.64 


1.35 


7078 


310 


35 


6.136 


0.75 


5.00 


3.229 


1.33 


3.06 


3.813 


0.36 


1.41 


1.981 


1.37 


3.25 



Notes. For each cluster, named in column (1), we provide the number of points in the surface brightness (2) and in the velocity dispersion (3) 
profile. For King and models, we list: the reduced best-fit photometric chi-squared;^ (Col. (4) and (10)), the mean (Col. (5) and (11)) and 
maximum (Col. (6) and (12)) photometric residuals, the reduced best-fit kinematic chi-squared Jj; (Col. (7) and (13)), the mean (Col. (8) and (14)) 
and maximum (Col. (9) and (15)) kinematic residuals. 



Table 5. Two-sided confidence intervals for the reduced ^-distribution with n degrees of freedom. 



Photometric Fits Kinematic Fits 



NGC 


n 


King^ 






^sup 


n 


King^ 






— ~2 

rC k.,Slip 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


104 


229 


3.49 


6.43 


0.85 


1.16 


15 


7.41 


10.37 


0.48 


1.67 


288 


83 


1.25 


3.89 


0.76 


1.27 


5 


0.44 


2.04 


0.23 


2.21 


362 


237 


3.11 


1.56 


0.85 


1.16 


7 


1.31 


3.35 


0.31 


2.01 


2419 


135 


1.98 


1.49 


0.81 


1.21 


5 


1.34 


0.47 


0.23 


2.21 


3201 


78 


1.31 


1.29 


0.75 


1.28 


15 


1.78 


4.01 


0.48 


1.67 


5139 


70 


3.75 


21.74 


0.74 


1.29 


36 


1.97 


3.41 


0.65 


1.42 


6121 


226 


1.46 


1.71 


0.85 


1.16 


9 


0.45 


0.58 


0.37 


1.88 


6218 


141 


1.19 


2.66 


0.81 


1.20 


10 


0.58 


0.77 


0.39 


1.83 


6254 


160 


5.05 


4.37 


0.82 


1.19 


5 


0.61 


1.84 


0.23 


2.21 


6341 


116 


8.44 


20.59 


0.79 


1.23 


7 


1.42 


2.35 


0.31 


2.01 


6656 


141 


1.02 


1.06 


0.81 


1.20 


6 


0.94 


1.70 


0.27 


2.10 


6809 


112 


1.17 


4.40 


0.79 


1.23 


12 


1.10 


2.97 


0.44 


1.75 


7078 


308 


6.14 


3.81 


0.87 


1.14 


34 


3.23 


1.98 


0.64 


1.43 



Notes. For each cluster, named in column (1), separately for the photometric and the kinematic fits, we provide the number of degrees of freedom 
of each fit (Col. (2) and (7)), the reduced best-fit chi-squared for King (Col. (3) and (8)) and for / <v) (Col. (4) and (9)) models, and the lower 
(Col. (5) and (10)) and upper (Col. (6) and (11)) boundaries of the two-sided 90% confidence level interval for the reduced ^-distribution with n 
degrees of freedom. 



Within the class of relaxed globular clusters, for NGC 362 
and NGC 7078 it is evident that the King models cannot repro- 
duce the observed surface brightness profiles, and that the f M 
models perform better, especially for describing the outer parts 
of the cluster. As to the observed velocity dispersion profiles, 
we see that for NGC 7078 the profile is formally more ad- 
equate (at the 99.9% confidence level, which, in the following, 
we denote by CL), while for NGC 362 the King profile is the 
closer to the observations (at the 90% CL). We should recall that 
these two globular clusters are flagged in the Harris catalogue 
as post-core-collapse clusters, although the value of the concen- 
tration parameter C is smaller than 2.5. Indeed, the observations 
indicate that there are some processes that cannot be captured 
by King models (with particular reference to the shallow cusp in 
the photometric profile of NGC 7078). By looking at the plots 
in Fig. CD NGC 104 and NGC 6341 appear to have both ob- 
served surface brightness and velocity dispersion profiles well 
represented by King models (for NGC 6341 the King falls 



within the 90% CI), in spite of the fact that, from Table[5] the fits 
cannot be considered entirely satisfactory. Curiously, NGC 6121 
is equally well described by the two models; the values of the 
reduced chi-squared are slightly lower for the King model, but 
in practice the quality of the two fits is similar. We notice that 
for NGC 104 and NGC 6121 it is particularly evident that the 
last point in the velocity dispersion profile is significantly higher 
than expected by both models. A partial explanation of this fact 
could be the very large extent of the radial interval in which it 
is calculated. Therefore, it is important to obtain more velocity 
data-points in the outer regions to clarify this issue. 

Best-fit profiles for globular clusters in the intermediate re- 
laxation class are shown in Fig. For NGC 288, NGC 6218, 
and NGC 6809 King models provide a better fit to both the sur- 
face brightness and the velocity dispersion profile (at the 90% 
CL in all cases). For the other globular clusters in this relaxation 
class the results are less sharp. In fact, for NGC 3201 and NGC 
6254 the surface brightness profiles are formally better repro- 
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NGC 362 - Surface Brightness Profile 



NGC 362 - Velocity Dispersion Profile 
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NGC 104 - Surface Brightness Profile 
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NGC 104 - Velocity Dispersion Profile 
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NGC 6121 - Surface Brightness Profile 
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NGC 6121 - Velocity Dispersion Profile 




NGC 6341 - Surface Brightness Profile 
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NGC 6341 - Velocity Dispersion Profile 
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Fig. 1. Fits by King models and anisotropic / (v) models to the surface brightness profiles and to the line-of-sight velocity dispersion 
profiles of relaxed globular clusters. In all panels, solid lines correspond to the King-model fits, dotted lines to / (v) -model fits; the 
vertical solid line marks the position of the King model projected half-light radius , R e , the dotted on e the position of the f (y} model 
projected half-light radius, R e . For the surface brightness profiles, the data from lTrager et al.l (1 19951) are indicated with circles, the 
data from other sources (see Sect. 13. U with squares. For each data-point, errors are shown as vertical error bars; in the case of the 
velocity dispersion profile, the horizontal bars indicate the length of the radial bin in which the data-points have been calculated and 
have no role in the fitting procedure (see AppendixlBt. 
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NGC 6218 - Surface Brightness Profile 
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NGC 6218 - Velocity Dispersion Profile 
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NGC 6656 - Velocity Dispersion Profile 
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NGC 3201 - Surface Brightness Profile 
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NGC 3201 - Velocity Dispersion Profile 
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NGC 6809 - Surface Brightness Profile 
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NGC 6809 - Velocity Dispersion Profile 
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NGC 288 - Velocity Dispersion Profile 
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Fig. 2. Fits by King models and anisotropic p** models to the surface brightness profiles and to the line-of-sight velocity dispersion 
grofiles of globular clusters in the intermediate relaxation condition, in the same format as in Fig.Q] 
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NGC 5139 - Surface Brightness Profile 



NGC 5139 - Velocity Dispersion Profile 
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Fig. 3. Fits by King models and anisotropic models to the surface brightness profiles and to the line-of-sight velocity dispersion 
profiles of partially relaxed globular clusters, in the same format as in Fig.Q] 



Table 6. Derived structural properties. 











King Models 














f' r} Models 








NGC 


c 


R c 




'•ti- 


M 


M/L 


Po 


Rc 


r a /r M 


m 


r a 


M 


M/L 


Po 


(1) 


(2) 


(3) 


(4) 


CS) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


104 


2.00 


22.60 


2.57 


2335.74 


7.181 


1.34 


5.011 


23.25 


1.783 


4.82 


8.59 


8.047 


1.50 


5.094 


288 


0.99 


79.49 


7.51 


896.85 


0.740 


1.88 


2.043 


58.24 


0.866 


7.53 


6.53 


0.850 


2.18 


2.427 


362 


1.77 


9.88 


2.65 


605.79 


1.867 


1.05 


4.826 


11.34 


1.786 


2.13 


3.80 


1.828 


1.03 


4.761 


2419 


1.42 


18.58 


19.56 


516.97 


7.843 


1.72 


1.874 


20.01 


0.958 


23.97 


22.95 


10.912 


2.40 


2.081 


3201 


1.30 


71.58 


3.85 


1532.62 


1.131 


1.91 


3.012 


72.31 


0.915 


5.14 


4.70 


1.088 


1.86 


3.026 


5139 


1.32 


127.68 


7.51 


2861.08 


26.446 


1.93 


3.537 


98.89 


0.978 


7.46 


7.30 


35.427 


2.58 


4.049 


6121 


1.62 


71.31 


2.62 


3144.17 


0.654 


1.10 


3.656 


72.45 


1.868 


4.40 


8.23 


0.750 


1.26 


3.661 


6218 


1.28 


47.81 


2.47 


982.45 


0.614 


1.96 


3.306 


43.45 


0.890 


3.03 


2.69 


0.786 


2.52 


3.577 


6254 


1.32 


50.02 


2.48 


1125.66 


1.532 


1.61 


3.741 


47.97 


0.564 


3.35 


1.89 


2.161 


2.67 


4.003 


6341 


1.69 


14.18 


2.12 


724.28 


2.866 


1.83 


4.638 


17.12 


1.523 


2.43 


3.70 


3.956 


2.53 


4.691 


6656 


1.38 


80.92 


3.13 


2057.81 


2.081 


1.11 


3.636 


82.71 


1.523 


4.53 


6.89 


2.337 


1.24 


3.652 


6809 


0.92 


110.09 


5.90 


1072.17 


0.604 


1.12 


2.214 


74.28 


0.867 


5.83 


5.06 


0.627 


1.16 


2.628 


7078 


1.86 


7.51 


1.70 


560.55 


3.976 


1.12 


5.207 


6.27 


1.793 


2.92 


5.24 


4.056 


1.14 


5.420 



Notes. For each cluster, listed in column (1), for the King models, we provide (2) the concentration index c = log(r tl /rn) (see Eqs. l lA.2t . JA.3H and 
(5) the truncation radius r tr , in arcsec; for the / <y) models, in Col. (10) the ratio between the anisotropy and the half-mass radius and in Col. (12) 
the anisotropy radius r a , defined as a(r a ) = 1 (see Eq. ( IB. 161 ) for the definition of a), in pc. For both King and f iv} models, we list: the core radius 
R c (defined in the standard way) in arcsec (Col. (3) and (9)), the intrinsic half-mass radius r M , in pc (Col. (4) and (11), respectively), the total mass 
M of the cluster (Col. (6) and (13)) expressed in units of 10 s Mq, the V band mass-to-light ratio in solar units (Col. (7) and (14)) and the logarithm 
of the central mass density p in Mq pc~ 3 (Col. (8) and (15)). 



duced by / (v) models (for the first, at 95% CL), while the cor- 
responding velocity dispersion profiles are formally better de- 
scribed by King models (at the 95% and 90% CL, respectively). 
For NGC 6656 both the surface brightness and the velocity dis- 
persion profile are approximately equally well reproduced by the 
two families of models (at the 90% CL in all cases). For NGC 
320 1 , NGC 6656, and NGC 6809 the velocity dispersion profiles 
have an irregular shape in the central regions: even if King mod- 
els formally perform better than / (v) models, they are unable to 
reproduce the observations. We tried to choose a different bin- 
ning for the data and we found that this irregularity does not 



depend on the way in which the observed velocity dispersion 
profile is constructed from the available data-set. Clearly, more 
data are necessary in order to obtain a more convincing descrip- 
tion of these systems. 

In Fig. [3] we show the best-fit profiles for partially relaxed 
globular clusters. For these two clusters we see discordant re- 
sults: for NGC 2419 / <v) models are more adequate for describ- 
ing the data (at the 99.99% CL and 90% CL for the photometric 
and kinematic fit, respectively), while for NGC 5139 King mod- 
els provide a better fit to the observed profiles. Even if formally 
King models perform better in describing the kinematic profiles 
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Fig. 4. Comparison between the values of the structural parameters derived from King and / (v) best-fit models. The top left panel 
shows, as crosses, the values of half-mass radii, the top right panel the total masses, the bottom left panel the central mass density 
and the bottom right panel the mass-to-light ratios. In each panel, the x-axis refers to the values calculated with King models, and 
the y-axis to those calculated with / (v) models; the dotted diagonal line denotes the identity relation. The intrinsic quantity ru would 
be best given in parsecs, but it is reported here in arcseconds for easier comparison with the observed profiles and with the projected 
half-light radius R e represented in the previous figures. 



of NGC 5139, they do not provide a satisfactory description of 
the kinematics of the central parts of the cluster (see Fig. [3); in 
this respect, the / (v) models give a better representation of the 
inner kinematics. 

To summarize, we found that, as expected, / (v) models tend 
to perform globally better than King models for the least relaxed 
globular cluster of our sample, NGC 2419. For NGC 2419, the 
good performance of / (v) models might correspond to the par- 
tial relaxation condition of the cluster, consistent with the phys- 
ical picture that motivates the definition of the / (v) models, as 
outlined in the Introduction. In addition, for three clusters in the 
second relaxation class (NGC 3201, NGC 6254,and NGC 6656), 
/ (y) models are competitive with King models. Furthermore, / (v) 
models can describe well a relatively steep central slope of the 
velocity dispersion profile, even when the corresponding photo- 
metric profile is cored, while King models and other isotropic 
truncated models (such as Wilson models) are unable to repro- 
duce this kinematical behavior. This fact is evident from the 
kinematic fits for NGC 2419, NGC 5139 and for possibly one 
cluster in the intermediate relaxation condition (NGC 6218). 

As far as the behavior of the photometric profiles at large 
radii is concerned, we see that, especially for NGC 104, 



NGC 362, NGC 6254, NGC 6341, and NGC 7078 King 
models do not provide a good description of the trunca- 
tio n as noted in a number of previous studies (in particular , 
see iMcLaughlin & van der Marell 120051: Uordi & Grebel 120101: 
Kflp per et al.1 l2010l in which the outermost parts of the sur- 
face brightness profiles are appropriately modeled by using N- 
body simulations). In some cases, the observed profile falls be- 
tween the King and the f (y} profiles. This suggests that trun- 
cated / (v) models might behave systematically better than King 
models for describing these stellar systems. To a large extent, 
the modification by truncation in phase space of the anisotropic 
(nontruncated) / (v) models is complementary to the generaliza- 
tion of the isotropic (truncated) King models to models char- 
acterized by anisotropic pressure, that is, the so-called Michie- 
King models (in which the truncated Maxwellian is asso ciated 
with the anisotropic factor o f the Eddington models; see[Michie 
119631 iGunn & GriffinlfT979l) . At this stage, it would be inter- 
esting to compare the behavior of the outer photometric and 
kinematic profiles of the two families of anisotropic truncated 
models, to evaluate the interplay between truncation and pres- 
sure anisotropy in the two different cases. Of course, the simple 
physical picture offered by King, / (v) , and Michie-King mod- 
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els still suffers from a number of limitations, as discussed in the 
Introduction. Such a simple picture is bound to fail in the model- 
ing of clusters in which core collapse has taken place. Therefore, 
it is not surprising that clusters such as NGC 362 and NGC 7078, 
even if they belong to the first relaxation class, are not described 
by isotropic King models as well as expected. 

In closing, we wish to reiterate that, in general, the kinematic 
fits are crucial to assess if a model is actually able to describe a 
given globular cluster. Unfortunately, the observed velocity dis- 
persion profiles are generally less accurate and less reliable, with 
respect to the surface brightness profiles; not only the outer parts 
(radii close to the truncation radius), but also the inner parts (in- 
side the half-light radius) are often not well sampled as would 
be desired. The present study confirms that in the future it would 
be desirable to acquire new and better kinematic data. 

4.2. King models vs. / (v) models 

The values of the relevant structural parameters derived from 
the best-fit models are presented in Table [6] To compare quanti- 
tatively the properties of the best-fit models selected in the two 
families of models for each globular cluster, we can correlate the 
values of the derived parameters, such as the half-mass radius 
tm, the total mass M, the central mass density po and the mass- 
to-light ratio M/L. These correlations are illustrated in Fig. [4] 

For the majority of the globular clusters considered in our 
sample, the values of the half-mass radius from f iy) models are 
larger than those obtained from King models; only for NGC 
362 the opposite is true. For NGC 288, NGC 5139 and NGC 
6809 the values calculated from the two models basically co- 
incide. By comparing the values of the total mass and of the 
central mass density calculated with the two families of mod- 
els, we see that for the whole sample (with few exceptions) the 
values calculated with / (v) models are larger than those calcu- 
lated with King models; this fact is not at all surprising, since 
the f iv> models are not truncated. As to the mass-to-light ratios, 
we see that there is not a tight correlation between the values 
calculated with King and / (v) models, the latter being almost al- 
ways larger. Similar trends are noted also in the structural prop- 
erties derived from (isotropic) models characterized by a more 
spatially extended truncation, such as the Wilson models (see 
iMcLaughlin & van der Marelll2005h . 

In Table [7] we list the values of the core and half-mass re- 
laxation times calculated by using the two best-fit dynamical 
models for each globular cluster; the clusters are listed in order 
of increasing King core relaxation times, and the separation in 
three classes of relaxation here adopted is marked with horizon- 
tal lines. Figure [5] illustrates the correlation between the values 
of these parameters. 

When considering the core relaxatio n times (calculated ac- 
cording to Eq. (10) of Djorgoyskl dl993h . as in the Harris 2010 
catalog), we see that the original division in the relaxation 
classes proposed in this paper on the basis of the values listed 
in the Harris catalog is confirmed. The only exception is NGC 
6341: according to the value of the core relaxation time calcu- 
lated with its best-fit / (v) model, this cluster should belong to 
the second class, rather than to the first. When we list globular 
clusters according to increasing core relaxation times estimated 
with the models identified in this paper, we see that the order of 
relaxed globular clusters changes with respect to that of Table Q] 
However, there is a general agreement between values of these 
quantities calculated with the two families of models. 

To calculate the half- mass relaxation times we followed the 
definition of Eq. (5) in Soi tzer & Hartl d 197 lh . which is based 



Table 7. Core and half-mass relaxation times for the best-fit 
King and p v ' models. 





King 


Models 




Models 


NGC 


logT c 


logr M 


io g r c 


logr M 


(1) 


(2) 


(3) 


(4) 


(5) 


362 


7.565 


8.971 


7.713 


8.825 


7078 


7.620 


8.820 


7.491 


9.176 


6121 


7.818 


8.774 


7.836 


9.137 


104 


7.846 


9.198 


7.922 


9.629 


6341 


7.881 


8.904 


8.142 


9.052 


6218 


8.140 


8.725 


8.142 


8.902 


6254 


8.269 


8.892 


8.334 


9.150 


6656 


8.418 


9.099 


8.451 


9.361 


3201 


8.523 


9.123 


8.545 


9.305 


6809 


8.835 


9.289 


8.529 


9.288 


288 


8.969 


9.482 


8.750 


9.509 


5139 


9.515 


10.140 


9.429 


10.191 


2419 


9.812 


10.537 


9.970 


10.731 



on the half-mass radius. We see that the values from p v > mod- 
els are larger than those from King models, except for NGC 362 
and NGC 6809 (for the latter cluster, the estimated values are 
approximately equal). In contrast with the case of the core re- 
laxation times, the values of the half-mass relaxation time turn 
out to be more model-dependent. Such dependence is likely to 
be due to the differences between the density profiles of the 
two families of models, which, in general, are less significant 
in the central regions and become more evident at radii larger 
than the half-mass radius. The introduction of a truncation for 
/ (v) models would also lead to different values of these pa- 
rameters. We notice that usually the half-mass relaxation time 
is calculated by inserting in the relevant definit ion directl y the 
(projected) half-light radius (see, f or example, iHarrisluOlOl and 
McLaughl in & van der Mare]|2005l) . By following the same pro- 
cedure, we found that the values of the half-mass relaxation 
times are less model dependent, and closer to the values listed 
in the Harris catalog. 

4.3. Comparison with previous studies 

In the case of King models, it is interesting to compare the val- 
ues of the structural parameters found in the present paper with 
those obtained in previous studies, as summarized in Table [8] 
By combining the formal errors from both analyses, th e com- 
parison of this paper with McLaughli n & van der MareH |2005) 
is the most significant because their values result from a fitting 
procedure which is similar to the one we have followed; for this 
reason, it is also, to some extent, the most surprising. We recall 
that in t his paper the errors on the paramete rs can be found in 
Table l3l iMcLaughlin & van der Mare! (120051) list the formal er- 
rors on the parameters in their Tables 10 and 12. No similarly 
detailed comparison with other papers could be made, because 
in general error analysis is not provided. 

We notice that for ten globular clusters (all but NGC 2419, 
NGC 6254, and NGC 70 78) our values of ¥ agree, wi t hin th e 
errors, with the values of McLaughlin & van der Marel 12001 : 
for the scale radius, agreement is found for eight objects (all 
but NGC 2419, NGC 6121, NGC 6254, NGC 6341, and NGC 
7078). We may argue that the discrepancy is partly due to the 
fact that our surface brightness profiles usually contain a larger 
number of data-points with respect to those of the cited paper 
(se e Sect. 13. Ill, even w hen we do not merge the original profiles 
by [Trager et alj d 19951) with those of other sources. 
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Fig. 5. Comparison between the values of the relaxation times derived from King and f (v) best-fit models. The left panel shows the 
values of core relaxation times, the right panel the half-mass relaxation times. The format is the same as in Fig. [4] 



Table 8. Comparison between the values of structural parameters from best-fit King models found in this paper and in previous 
studies. See Table[3]for the description of column entries. 



NGC 


Ref. 


*P 


ro 


A<o 


V 


104 


(0) 


8.58 


23.09 


14.33 


12.27 




(1) 


9.33 


15.25 


14.42 


1 1.5 




(2) 


8.70 


22.84 


14.42 






(3) 


8.6 


23.19 


14.36 


15.27 




(4) 


8.81 


22.01 


14.38 


11.0 


362 


(0) 


7.80 


10.21 


14.66 


8.31 




(1) 


8.93 


7.19 


14.79 


6.4 




(2) 


6.80 


10.79 


14.79 






(3) 


7.9 


10.05 


14.66 


11.12 




(4) 


7.76 


11.16 


14.80 


6.4 


3201 


(0) 


6.17 


76.99 


18.35 


4.28 




(1) 


6.17 


86.22 


18.96 


5.2 




(2) 


6.21 


93.55 


18.96 






(3) 


6.1 


77.68 


18.30 


4.21 




(4) 


6.14 


83.98 


19.00 


5.0 


6121 


(0) 


7.32 


74.56 


17.00 


4.01 




(1) 


7.24 


45.38 


17.88 


4.2 




(2) 


7.20 


52.37 


17.88 






(3) 


7.4 


72.57 


16.94 


5.25 




(4) 


7.40 


72.42 


17.95 


4.0 


6254 


(0) 


6.26 


53.63 


16.88 


6.21 




(1) 


6.55 


48.45 


17.70 


6.6 




(2) 


6.55 


54.49 


17.70 






(3) 


6.5 


49.41 


16.85 


6.17 




(4) 


6.48 


49.19 


17.70 


6.6 


6656 


(0) 


6.47 


86.18 


16.41 


6.47 




(1) 


6.17 


83.19 


17.40 


9.0 




(2) 


6.21 


91.41 


17.40 






(3) 


6.5 


85.15 


16.38 


8.56 




(4) 


6.48 


84.96 


17.42 


7.8 


7078 


(0) 


8.09 


7.72 


14.07 


11.83 




(1) 


12.32 


1.96 


14.21 


12.0 




(2) 


10.75 


4.09 


14.21 






(4) 


9.72 


8.49 


14.21 


13.5 



NGC 


Ref. 




ro 


y"0 


V 


288 


(0) 


4.82 


91.03 


20.02 


2.85 




(1) 


4.33 


94.05 


20.00 


2.9 




(2) 


4.65 


98.41 


20.00 






(3) 


4.8 


92.87 


20.00 


2.79 




(4) 


4.80 


92.79 


20.05 


2.9 


2419 


(0) 


6.62 


19.70 


19.43 


5.04 




(1) 


6.55 


19.88 


19.77 


3.0 




(2) 


6.55 


22.19 


19.77 






(3) 


6.5 


20.60 


19.44 


5.32 




(4) 


6.44 


20.47 


19.67 


4.0 


5139 


(0) 


6.27 


136.94 


16.42 


14.83 




(1) 


5.77 


157.85 


16.81 


16.0 




(2) 


5.94 


168.01 


16.81 






(3) 


6.2 


141.20 


16.44 


13.75 




(4) 


6.21 


152.73 


16.81 


16.8 


6218 


(0) 


6.11 


51.56 


17.65 


3.93 




(1) 


6.55 


37.57 


18.00 


4.5 




(2) 


6.48 


42.38 


18.00 






(3) 


6.1 


51.74 


17.62 


4.66 




(4) 


6.33 


50.70 


18.10 


4.5 


6341 


(0) 


7.54 


14.72 


15.31 


9.28 




(1) 


8.22 


11.27 


15.46 


5.9 




(2) 


7.92 


14.56 


15.46 






(3) 


7.5 


16.15 


15.61 


8.71 




(4) 


7.50 


16.20 


15.47 


6.0 


6809 


(0) 


4.44 


129.11 


19.12 


2.92 




(1) 


3.17 


192.94 


19.40 


4.9 




(2) 


3.53 


215.54 


19.40 






(3) 


4.5 


126.41 


19.07 


3.73 




(4) 


4.49 


126.22 


19.36 


4.0 



Notes. For NGC 7078 the entry marked with (3) is missing: the authors of this work decided to exclude this globular cluster from their analysis 
because it is core-collapsed. In references (2) and (4) only the core radii are given; therefore we used the values of R c and T to calculate ro (see 
App.rAj}. 

References. (0) This work; (D lPrvor & Mevlanll 19931 ; (2) lTrager et aflll995l ; (3) lMcLaughlin & van der Marelll2005l ; (4) lHarrisll201d . 
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The discrepancies between the values of the central surface 
brightness derived in this paper and those resulting from pre- 
vious studies are primarily due the correction for the extinc- 
tion, usually neglected, introduced in our analysis, following 
iMcLa ughlin & van der M arel (2005). In fact, such discrepancies 
are not severe; only three objects (NGC 6341, NGC 6809, and 
NGC 7078) have values which are not consistent, within the er- 
rors, with respect to values determined in previous studies. The 
reason for this is the fact that we ad ded to these profiles the mor e 
recent and more accurate data from|Noyola & Gebhardjj d2006h . 

The structural parameter for which the differences between 
the values obtained in our analysis and those in the literature are 
the most relevant is the central line-of-sight velocity dispersion. 
For NGC 288, NGC 2419, NGC 3201, NGC 6121, NGC 6254, 
and NGC 7078 we see that our values agrees, within the errors, 
with at least one of the values found in the litera ture (only four 
of them with lMcLaughlin & van derM arel 2005). This discrep- 
ancy is not at all surprising, because the kinematic data are the 
most uncertain. This fact gives one further argument for the need 
for more numerous and more accurate kinematic data for these 
systems. 

4.4. Central slope of the photometric and kinematic profiles 

When considering the surface brightness profiles of our selected 
globular clusters, we note that there are cases in which the ob- 
served photometric profiles deviate from the calculated ones at 
small radii. In particular, by focusing on the innermost regions 
of the profiles, we see that this is the case for NGC 6121, NGC 
6218, and NGC 7078, for which the models are underluminous, 
and NGC 288, NGC 6656, and NGC 6809, in which the models 
are overluminous, and NGC 3201, for which the observed cen- 
tral surface brightness appears to oscillate; in spite of these local 
discrepancies, the global values of the statistical indicators may 
be satisfactory (see Table|5]l. 

As far as the velocity dispersion profiles are concerned, in 
four cases (NGC 288, NGC 3201, NGC 6121, and NGC 6656), 
both models overpredict the central data-points, while, as men- 
tioned in Sect. ED four clusters (NGC 2419, NGC 5139, and, 
possibly, NGC 6218 and NGC 6254) show a relatively large gra- 
dient of the profile in the central regions. At variance with the 
study of elliptical galaxies, the kinematic profiles of globular 
clusters are often undersampled inside the half-light radius. 

Recently, the central cusps in the observed photometric and 
kinematic profiles of some globular clusters have been inter- 
preted as clues of the presence of an intermediate-mass black 
hole (IM BH) in the center of the system (see the analytical 
model by Bahcall & Wolf 1976 a nd the N-body s imulations by 
Baumgar dt et al.l 120051 andlNovola & Baumgardtl l201 ll) . and a 
variety of dynamical models (either defined from distribution 
functions or as solutions of the Jeans equations), have been 
used in order to constrain the mass of such central object (the 
best k nown example is th e c ontroversial case of u> Cen, stud - 
ied by iNovola et afll2008l and Ivan der Marel & Andersoiill20T 
with different conclusions). However, [Vesperini & Trenti 1'. 



Table 9. Global anisotropy parameter for the best-fit models. 



, SLUU- 

2010, 

:20i0) 



showed that these shallow photometric cusps are not decisive 
signatures of the presence of an IMBH, and that they might be 
related to other dynamical processes; moreover, the authors em- 
phasize the fact that the typical accuracy in the data may be in- 
sufficient to characterize the slope of the profile as desired. 

As shown in the previous sections, the presence of radially- 
biased pressure anisotropy (which occurs in the contex t of the 
famil y of / (v) models or of the Michie-King models, see lMichia 
119631) can also produce a relatively rapid decline in the central 



Relaxed 


Intermediate 


Partially Relaxed 


NGC 


K 


NGC 


K 


NGC k 


362 


1.353 


6218 


1.835 


5139 1.754 


7078 


1.301 


6254 


2.307 


2419 1.771 


104 


1.301 


6656 


1.443 




6121 


1.320 


3201 


1.810 




6341 


1.443 


6809 


1.859 








288 


1.859 





part of the velocity dispersion profile. Therefore, here we are 
reiterating a point already noted in the literature, that the pres- 
ence of a central IMBH should not be considered as the only 
physical explanation of the existence of a central kinematical 
peak0. Unfortunately, our independent conclusion only confirms 
that the interpretation of this interesting kinematical feature is 
more model dependent than desired. We recall that / (v) models 
ar e characterized by a " realistic" anisotropy profile (see Fig. 6 
in iTrenti & Bertir] [2005): the central regions are more isotropic 
than the outer ones in velocity space, because the models rep- 
resent a scenario in which violent relaxation has acted more 
efficiently in the center. This differential kinematical feature is 
not always present in dynamical models based on the Jeans ap- 
proach, which, to obtain a fast decline in the velocity dispersion 
profile in the absence of a central IMBH, usually requires very 
high values of the anisotropy parameter, even in t he central re- 
gions of the cluste r (see, for example, Sect 4.2 in INovola et al.l 
120081 or Sect 5.3 in Lutzgendorf^talJ|2011|). 



4.5. Radial orbit instability 

Systems in which there is a great amount of radial kinetic energy 
with respect to tangential kinetic energy are subject to the radial 

orb it instability. 

iPolvachenko & Shukhmanl (Il98ll) (see also 
Fridman & P olvachenkol 1 19841) introduced a parameter de- 
fined as: 



2K r 

Kj 



(1) 



where K r and Kj represent the radial and the tangential com- 
ponent of the total kinetic energy, respectively, and argued that 
when: 



k> 1.7 + 0.25 



(2) 



radial orbit instability occurs. Actually, different families of 
models generally have different threshold values for the instabil- 
ity. A detailed discussion of the onset of the radial orb it instabil- 
ity for t he family of / (v) m odels ca n be found in | Trenti & Bertinl 
(120051) . ITrenti et al.l (120051) . and ITrenti & Bertinl d2006l) . where 
the validity of the criterion expressed in Eq. (|2]i is discussed. 

We calculated the values of the parameter k for the / (v) mod- 
els that provide the best fit to the observed profiles of the glob- 
ular clusters of the sample (see Table|9]l. It is interesting to note 
that globular clusters belonging to a given relaxation class tend 
to have similar values of the stability parameter (k a 1.3, 1.8, 1.7 
for the first, second, and third relaxation class, respectively; in 



5 Of course, / (v) models are able to reproduce a slope in both the 
photometric and kinematic profile, only when high values of the con- 
centration parameter are considered. 
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this respect, NGC 6656 and NGC 6254 appear to exhibit an ex- 
ceptional behavior). In other words, the relaxed class is found to 
be more isotropic by the f v) diagnostics. 

The majority of the globular clusters that have k ^ 1.7 are 
likely to be in a condition of marginal instability. The case of 
NGC 6254 does remain problematic. 

We argue that the introduction of a truncation in phase space 
to the family of / (v) models might have a stabilizing effect, since 
such truncation will affect primarily the outer parts of a given 
configuration, which are dominated by radially-biased pressure 
anisotropy. Therefore the truncation is likely to reduce the global 
value of the radial component of the total kinetic energy. 

This is one more reason to undertake the study of truncated 
/ (v) models for application to globular clusters. We plan it for 
the immediate future. 

5. Conclusions 

In this paper we have performed a detailed combined photomet- 
ric and kinematic study of a sample of Galactic globular clusters, 
representing systems under different relaxation conditions. For 
these objects, surface brightness and velocity dispersion profiles 
have been fitted by means of two different families of dynam- 
ical models, the truncated isotropic King models and the non- 
truncated, anisotropic / <v) models. The analysis has been carried 
out by following the same procedure used in the past to study 
the dynamics of elliptical galaxies. Each globular cluster is then 
associated with two best-fit models. The main conclusions can 
be summarized as follows: 

- The expected trend, that King models should perform better 
for more relaxed globular clusters, has been checked to exist 
but it is not as sharp as anticipated. The two clusters (NGC 
104 and NGC 6341) for which the global fit by King mod- 
els is most convincing indeed belong to the class of relaxed 
objects. King models tend to offer a good representation of 
the observed photometric profiles (as is commonly reported), 
regardless of the relaxation condition of the system (but a 
statement of this kind should also be supported by the rele- 
vant statistical indicators; see Table|5]l. However, the quality 
of the fits by King models to the kinematic profiles remains 
to be proved, even for relaxed clusters, because of the few 
data-points and the large error bars in the observed profiles. 
Three clusters for which the King models appear to be inad- 
equate do not actually come as a surprise: NGC 2419 is the 
least relaxed cluster of the sample and NGC 362 and NGC 
7078 are suspected to be post-core-collapse clusters. 

- The second expected trend, that less relaxed clusters might 
exhibit the characteristic signature of incomplete violent re- 
laxation, is also partly present but is not as sharp as might 
have been hoped for. Some cases indeed point to a signifi- 
cant role of radially-biased pressure anisotropy. The least re- 
laxed cluster, NGC 2419, is well described by the / (v) mod- 
els. For the second least relaxed cluster, NGC 5139, the cen- 
tral shallow cusp in the velocity dispersion profile appears 
to be well captured by the models. A marginal indica- 
tion in favor of the / (v) models also comes from inspection 
of the inner kinematic profiles of NGC 62 1 8 and NGC 6254, 
although these two clusters are not among the least relaxed 
objects. In contrast, King models and other isotropic mod- 
els (such as the spherical Wilson models) have difficulty in 
matching significant velocity gradients inside the half-light 
radius. Therefore, the partial success of the / (v) models sug- 
gests that for some globular clusters radially -biased pressure 



anisotropy may be important. This property could be exam- 
ined further by means of better spatially -resolved kinematic 
data in the inner regions (i.e., at radii out to approximately 
the half-light radius). This result is in line with the con clu- 
sions of recent papers: in particular, see Ibata et al. 20Q]and 
references therein, based on the application of King-Michie 
models. 

- In some clusters, regardless of the relaxation condition, some 
qualitative characteristics of the observed profiles are missed 
by both families of models considered in this paper. It may 
be that part of these cases would be resolved by a study in 
terms of truncated / (v) models. But it may also be that other 
ingredients ignored in the paper, such as rotation (solid-body 
or differential), play a role. 

- The paper demonstrates that the values of some structural 
parameters, such as the total mass and the half-mass radius, 
can be significantly model-dependent. In view of the results 
listed in the previous items, this is a clear warning against 
an indiscriminate use of structural parameters for globular 
clusters based on only one family of models (the spherical 
King models). 

- In general, the kinematic fits are crucial to assess if a model 
is actually suited to describe a given globular cluster. The 
main issue in testing dynamical models on globular clusters 
is therefore the general lack of good kinematic data: the data 
are available for a small fraction of the population of Galactic 
globular clusters and generally made of a small number of 
data-points, not well distributed in radius. Surprisingly, the 
kinematic profile is often not well sampled as desired inside 
the half-mass radius. As discussed in the paper this is a key 
region for confronting the performance of different dynami- 
cal models. In addition, accurate data in the outermost parts, 
close to the truncation radius,would touch on other important 
issues, such as the role of tides. Only after the acquisition of 
good kinematical profiles will it be possible to address prop- 
erly the issue of dark matter in globular clusters. 
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Appendix A: The dynamical models 

A.1. Spherical isotropic King models 

The lKingl(ll966l) one-component, spherical, isotropic models are 
based on the following distribution function: 



_ . A [exp(-a£') - exp(-a£Vj)] E < Eq 
/k ~ i E > E Q 



(A.l) 



where A, a, Eq are positive constants, defining two scales and 
one dimensionless parameter, and E represents the specific en- 
ergy E = v 1 12 + O(r), where <J>(r) is the mean-field gravitational 
potential, to be determined from the Poisson equation. The quan- 
tity £0 is a threshold energy, above which the stars are consid- 
ered unbound; it can be translated into a truncation radius, r tr , 
for the system. 

We recall that the radial scale ro can be defined as 



ro 



\A-nGapQ 



1/2 



where po is the central mass density. The other dimensional scale 
is the total mass of the cluster or, as an alternative, its central 
velocity dispersion. 

The dimensionless concentration parameter is given by the 
central dimensionless potential *P = a[0(r tr ) - ^(O)] or, alterna- 
tively, by the index 



(A.3) 



c=log U ' 



because these two quantities are in a one-to-one relation. 

The quantity listed in the Harris catalog as concentration pa- 
rameter, in this paper indicated with C, is defined using the stan- 
dard core radius R c in place of ro in Eq. ( IA.3I ). The values of C 
for a model identified by *F 4 are slightly larger than the corre- 
sponding values of c; in fact, for these models 0.8 ;S R c /ro 55 1. 

A.2. Anisotropic f M models 

Several families of dynamical models have been developed to 
represent the final state of numerical simulations of the violent 
relaxation process thought to be associated with the formation of 
bright elli ptical galaxies via collisionless collapse (for a review, 
see lBertin & Stiavellilll993l) . These models show a characteris- 
tic anisotropy profile, with an inner isotropic core and an outer 
envelope that becomes dominated by radially-biased anisotropic 
pressure. They provide a good representation of the photometric 
and kinematic properties of elliptical galaxies. Here we will refer 
to the family of spherical, anisotropic, non-trunc ated / (v - > models 
(whic h have been revisited recently in detail bv lBertin & Trentil 
120031) . 

The distribution function that defines these models depends 
on specific energy E and angular momentum J: 



/( v) = 



A exp 

v 



a E d ^ |£| 3/2 ) 



E < 
E > , 



(A.4) 



where A, a, d, and v are positive constants, defining two scales 
and two dimensionless paramet ers. For applications, as de- 
scribed by Ber tin & Trentil ((2003), the dimensionless parameter 
v can be fixed at v = 1 . Therefore, similarly to the King models, 
after integration of the relevant Poisson equation the / (v) mod- 
els are a one-parameter family of models, parametrized by their 
central concentration, which can be expressed by the central di- 
mensionless potential ¥ = -a<E>(0). 

The physical scales can be expressed as: 



r SC aie = <T 1/V a- I/4 
M scale = d-^a-V* 



(A.5) 



(A.2) 



By definition, these models are non-truncated and, because 
of this, are likely to be less suited to describe the outer parts of 
globular clusters. A study of globular clusters based on truncated 
/ (v ) models is postponed to a separate investigation. 

We recall that the surface brightness profiles for concentrated 
models (¥ ^ 7) are very close to de Vaucouleurs profile, while 
for low values of ¥ the models exhibit a sizeable core. 



Appendix B: Fitting procedure 

In the following we describe the procedure that we have adopted 
to perform the stati stical analysis of the data. We basically follow 
lBertinetal.1 (119881) . 
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B.1. Photometric fit 



B.2. Kinematic fit 



For each family of models considered in this paper, the photo- 
metric fit determines which equilibrium model has the projected 
mass distribution that best reproduces the surface brightness pro- 
file, under the assumption that: 



(B.l) 



where 'E(R) is the model surface mass density, A(R) the model 
surface luminosity density, and the mass-to-light ratio M/L is 
considered to be constant in the cluster. From the photometric 
data, we perform a fit that allows us to determine three parame- 
ters for each model: 

- W, the concentration parameter, which determines the shape 
of the surface brightness and velocity dispersion profiles; 

- r s , the scale radius, that is ro for the King models and r sca i e 
for the / (v) models; 

- fj.0, the central surface brightness. 

We might determine the values of these parameters by mini- 
mizing the quantity: 



A r m{Rd-S{R i lr % )+^' 1 



Z 



Sin/ 



(B.2) 



where S is the surface density expressed in magnitudes, but we 
found that it is more convenient to minimize: 



X 



l(Ri) - kZ(Ri/r s ) 



Sit 



(B.3) 



because the chi -squared function turns out to be mor e stable, as 
pointed out by iMcLaughlin & van der M arel The con- 

nection between Eqs. ( IB. 21 ) and ( IB.3b is easily made, by taking 
into account the following relations: 

m(Ri)-m = -2.51og[/(*/)] 

S(R i /r s ) = -2.5log\t(R i /r s )] (B.4) 
jiio - niQ = 2.5 log A: , 

where the zero-point mo = 26.422 allows us to express the ob- 
served surface brightness m{R[) in solar units; the observed lu- 
minosities l(Rj) are expressed in solar units and 61/ are calculated 
from 5nij and represent the errors on luminosities. The quantity 
2(/?,/r s ) is the surface density, normalized to its central value; 
from this point of view, the parameter k is related to the mass- 
to-light ratio. 

Even if we calculated the parameters from the luminosities, 
we report the results in terms of magnitudes, because in this way 
a comparison with the data is more natural. To calculate the sur- 
face brightness profile fi(R) that best reproduces the data, we 
choose the model identified by *P, we calculate the normalized 
projected mass-density profile, and then we rescale it radially 
with r s and vertically with 

Note that we have only two independent fit parameters, be- 
cause k = r s ). Hence, the value of the reduced^ is: 



N p -2 



(B.5) 



To have a quantitative estimate of the quality of the fit, it is con- 
venient to calculate also the mean and maximum residuals: 



1 1/2 



(A/i) n 



max i=1 ,...jv p \m(Ri) - fi(Ri)\ . 



The parameters identified by the photometric fit determine a 
model, characterized by the parameter *F, the profiles of which 
are rescaled with r s . At this point, we perform a fit to the kine- 
matic data to find the central line-of-sight velocity dispersion, 
that is the velocity scale needed to rescale vertically the nor- 
malized projected velocity dispersion profile calculated from the 
dynamical models. We calculate the value of this scale, V, as the 
one that minimizes: 



A'k 

z 



oiR,) - Vcr P (Ri) 



5cTi 



(B.7) 



where cr P (/?,) is the projected velocity dispersion calculated from 
the model, dependent on the rescaled radial coordinate. In this 
case, the values of the reduced x 1 an d of the residuals can be 
found from the following expressions: 



~2 



JVv-l 



2l 



1/2 



(B.8) 



(Air) = [ziiMRd - Vcr P (Ri)l 
(Acr) max = max^i,..,^ \<r(Ri) - Vcr P (Ri)\ . 

One might argue that a more sensible way to perform the fit 
would be by means of a single fit procedure, using the combined 
Xtot = xl + x\> defined as the sum of the photometric and of 
the kinematic contribution. We did perform tests of this different 
procedure and found that the results are equivalent to those ob- 
tained by performing the photometric fit first and then the kine- 
matic fit at fixed *T and r s . This confirms the qualitative expec- 
tation that the kinematical data, being less numerous and less 
accurate with respect to the photometric ones, have little weight 
in determining *P and r s and are only needed to determine the 
scale V, that is, the relevant mass-to-light ratio. 

Finally, we addressed the issue of whether a different rep- 
resentation of the results of the best-fit models would be more 
suitable than the one used in Figs.Q~|[3] Figure iBTI compares the 
appearance of a fit for a given object using a linear and a loga- 
rithmic representation of the radial scale. Based on this and other 
tests, we preferred to adopt a mixed representation, as used in the 
main text of the paper. 



B.3. Goodness of the photometric and kinematic fits 

To measure the goodness of the photometric and kinematic fits, 
we referred to the confidence intervals on the x 1 -distribution. 
The probability density function of the^ 2 distribution is defined 



f(x; n) = 



2"/ 2 F(n/2) 



(B.9) 



when x > 0, and zero otherwise. The parameter n corresponds 
to the number of degrees of freedom and T de notes the gamma 
function (e.g., see Abramowitz & Steguni n~972l Sect. 26.4). 

The two-sided confidence interval relative to a confidence 
level a on the y 2 probability density function is defined as the 
interval \x 2 mi ,x*] such that 



-t 



and that: 



f(x; n)dx ; 



/(*; n) dx = I f(x; n) dx 



(B.10) 



(B.ll) 
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B.4. Errors on the best-fit parameters 

For th e methods described below, we refer mainly to lPress et al.l 
(2007). In the case of the photometric fit, we used the following 
procedure to calculate the formal errors on the parameters. First 
we calculated the Hessian matrix: 



' ; dxjdxj 



(B.12) 



where i = 1,2,3 and x\ — *P, xi — r s and x^ = /xq; then we 
calculated the covariance matrix: 



[E] = 2[HY 



Finally, we obtained the errors: 6x[ = (En) 



(B.13) 



1/2 



In the case of the kinematic fit, in which depends analyti 
cally only on the parameter V, we can immediately calculate: 



6V = 



1/2 



(B.14) 



In addition to errors, we identified the relevant confidence 
regions and intervals for the various parameters, after defining a 
likelihood function A = e ~ x ~l 2 . First, we calculated confidence 



regions in the Q¥, r s ) plane, after marginalizing the likelihood 
function on k, in the following way: 



A(T,r,) 



(B.15) 



We identified the regions corresponding to confidence levels of 
68.3% and 95.4%. To check the results, we compared them with 
the regions defined by the curves of constant y 1 w ith the ap- 
propriate values (see Sect. 15.6 in iPress et alJl2007h . For each 
globular cluster and for the families of models, we notice that 
the regions identified with these two methods overlap in a con- 
sistent way. 

We also calculated the confidence intervals on the parameters 
*P, r s and V (for the last parameter, the calculation can be done 
immediately, because the kinematic likelihood depends only on 
V; for the other parameters, we have to marginalize once more 
the likelihood function), corresponding to confidence levels of 
68.3% and 95.4%. Comparing them with the errors calculated 
with the covariance matrix method, we found that the results are 
consistent. 

In Fig. IB.2l we show an example of the confidence regions 
and intervals calculated for the King best-fit model for the glob- 
ular cluster NGC 6121. The top left panel shows the overlap of 
the confidence regions (solid lines) corresponding to confidence 
levels of 68.3% and 95.4%, and the curves of constant^ 2 (dashed 
lines). It is clear that there is a good agreement between the two 
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Fig. B.2. Confidence regions and confidence intervals on NGC 6121 King model parameters. The top left figure shows the overlap of 
the confidence regions, corresponding to confidence levels of 68.3% and 95.4% (solid lines), and the curves of constant^- 2 (dashed 
lines); the black dot marks the position of the maximum of the likelihood (that is the minimum of x 2 )- The top right figure shows 
the confidence intervals for the likelihood depending on V; the bottom figures show the confidence intervals on the parameter *F 
and ro; the dark grey area corresponds to the confidence level of 68.3%, the light grey area to the confidence level of 95.4%; a solid 
vertical line marks the position of the maximum of the likelihood (that is the minimum of^ 2 ). 

relevant pairs of curves. The bottom panels show the confidence The ratio of this radius to the half-mass radius measures how 
intervals calculated on the marginalized likelihood for the pa- large is the radially-biased anisotropic part of the globular clus- 
rameters *P and ro, and the top right panel those for the parameter ter: the smaller r a , the larger the anisotropic region of the cluster. 
V; in these figures the confidence intervals are shown as shaded 
areas. By comparing the extent of the dark grey area with the 
values of the uncertainties on the corresponding parameters (see 
Table [3}, we see that the different methods lead to consistent re- 
sults. 

B.5. Derived parameters 

Once the best-fit parameters have been secured, we may calcu- 
late other quantities, which represent some important character- 
istics of globular clusters. 

In particular, for each family of models we calculated the 
core radius R c , that is the radial position where the surface 
brightness equals half its central value, the half-mass radius r^, 
the total cluster mass M, the central mass density po, the mass- 
to-light ratio Mj L, the core relaxation time log T c and the half- 
mass relaxation time, log 7m ■ In addition, for King models we 
calculated the concentration parameter c and the truncation ra- 
dius r tr . In turn, for the / (v) models we calculated the anisotropy 
radius r a , defined as the radial position where a(r a ) = 1, with 

a = 2 - «v 2 ) + <v 2 »/<v 2 > . (B.16) 
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